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SIGNAL-TO- NOISE  RATIOS  REQUIRED  FOR  SHORT-TERM 
NARROWBAND  DETECTION  OF  GAUSSIAN  PROCESSES 


INTRODUCTION 

Rules  of  thumb  are  available  for  approximating  the  required  signal -to-noise 
ratio  (SNR)  for  broadband  and  narrowband  energy  detection  of  signals  in  noise 
(see  reference  1).  However,  they  are  generally  applicable  or  appropriate  only 
when  certain  requirements  (long  observation  (averaging)  time,  for  example)  are 
satisfied.  Attempts  to  apply  these  rules  outside  their  often  ill-defined  ranges  of 
applicability  can  result  in  very  misleading  conclusions  about  system  perform¬ 
ance  and  capability.  This  is  particularly  true  for  high- resolution,  narrowband 
detection  systems  where  the  analysis  bandwidth  is  so  small  that  practical 
observation  times  lead  to  only  a  few  effectively  independent  samples  of  the 
processes  under  investigation;  this  is  called  short-term  detection  here. 

The  rules  of  thumb  may  be  inapplicable  because  they  are  often  based  on  a 
deflection  criterion  of  system  performance,  such  as  the  square  of  the  difference 
of  the  mean  outputs,  with  and  without  signal,  divided  by  some  variance  of  the 
system  output.  Since  these  rules  use  only  first  and  second  moments  of  the  sys¬ 
tem  output  decision  variable,  they  are  incomplete  statistical  descriptors. 
Furthermore,  the  two  lowest  moments  are  often  employed  in  a  Gaussian  approxi¬ 
mation  in  order  to  estimate  probabilities  of  detection  and  false  alarms.  This 
approach  can  lead  to  an  obvious  contradiction,  such  as  stating  that  no  value  of 
SNR,  no  matter  how  large,  will  yield  very  high  detection  probabilities  for  a 
limited  observation  time.  * 

The  best  way  to  avoid  such  situations,  and  the  approach  this  report  dis¬ 
cusses,  is  to  evaluate  the  exact  probability  distribution  of  the  decision  variable 
without  making  a  Gaussian  assumption  —  a  method  that  is  particularly  relevant 
in  short-term  detection  where  only  a  few  independent  samples  of  the  system 
output  are  available. 

*See,  for  example,  reference  1,  equations  (37)  through  (40),  where  deflection 
criteria  Pj  and  P2  can  not  exceed  TWB,  and,  in  fact,  the  required  S0/N0  —  » 
as  Plf  P2  -  TWS. 
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PROBLEM  DEFINITION 

The  detector  of  interest  here  is  depicted  in  figure  1.  The  input  x(t)  is 
composed  of  stationary,  zero-mean,  Gaussian  signal  s(t)  and  noise  n(t),  or 
noise  n(t)  alone.  It  is  assumed  that  the  narrowband  filter  (NBF)  is  centered 
on  the  signal  and  has  the  same  bandwidth,  B  Hz,  as  the  signal  spectrum.  (The 
effects  of  over-resolution  and  under-resolution,  and  of  mismatched  center  fre¬ 
quencies,  are  discussed  in  reference  1. )  The  input  noise  is  assumed  to  be 
fairly  flat  over  a  frequency  interval  wider  than  B  Hz.  The  sampler  takes  a 
sample  of  the  squared-envelope  (SE)S2(t)  approximately  every  B-1  sec,  * 
accumulating  these  SE  samples  for  an  observation  interval  of  T  sec.  The 
threshold,  A ,  with  which  z  is  compared,  is  fixed.  (The  effects  of  using  a 
tracking  threshold  based  on  a  few  noise-alone  frequency  bins  are  discussed  in 
reference  1  and  thoroughly  investigated  in  reference  2. ) 


SQUARE-LAW  ENVELOPE  DETECTOR 


Figure  1.  Narrowband  Detector 

If,  after  an  observation  interval  of  T  sec,  random  variable  (RV)  z  exceeds 
A  when  a  signal  is  absent,  a  false  alarm  occurs.  It  should  be  noted  that  this 
false  alarm  definition  applies  to  just  one  NBF,  not  to  a  bank  of  NBFs.  Also,  a 
false  alarm  can  occur  only  once  every  T  sec,  not  every  B-1  sec.  If  A  is  ex¬ 
ceeded  with  a  signal  present,  a  detection  occurs.  The  problem,  then,  is  to 
calculate  the  false  alarm  and  detection  probabilities  of  the  detector  in  figure  1 
as  a  function  of  the  input  SNR. 


♦This  is  a  reasonable  rule  for  analysis  filters  with  fairly  rectangular  pass- 
band  characteristics,  but  samples  should  be  taken  more  often  for  filters  with 
rounded  characteristics.  In  the  latter  case,  however,  the  SE  samples  are 
statistically  dependent.  (Appendix  A  provides  a  method  for  dealing  with  this 
dependency. ) 
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PROBLEM  SOLUTION 


The  process  y(t) 


aJ  the  NBF  output  has  power 


where 


9 


2  = 


N,  noise-alone  input  to  NBF 
S  +  N,  signal-plus-noise  input  to  NBF 


(1) 


That  is,  N  is  the  noise  power  at  the  NBF  output,  and  S  (if  present)  is  the  sig¬ 
nal  power  at  the  NBF  output.  These  are  the  powers  in  the  B  Hz  bandwidth  and 
are  not  spectral  levels  in  a  1-Hz  band. 


It  will  be  assumed  that  the  samples  every  B_1  sec  at  the  NBF  (or  detec¬ 
tor)  output  are  statistically  independent.  This  is  a  fair  approximation*  if  B 
is  interpreted  as  the  effective  bandwidth  of  the  signal  process  and  the  NBF; 
that  is, 


(2) 


where  Gy(f)  is  the  NBF  output  power  density  spectrum.  The  number  of  inde¬ 
pendent  SE  samples  in  observation  interval  T  sec  is,  then,  BT  +  1  and  will 
be  denoted  by  M.  (The  equivalent  number  of  independent  SE  samples  in  inter¬ 
val  T  for  statistically  dependent  SE  samples  is  considered  in  Appendix  A. ) 

The  generic  problem,  therefore,  is  that  of  calculating  the  cumulative  dis¬ 
tribution  of  the  RV  z, 

z  -  E  wk  *  (3) 

k  =  1 


*See  footnote  on  p.  2. 
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where  are  statistically  Independent  samples  of  the  envelope  of  a  narrow- 
band  Gaussian  process.  Random  variable  z  is  a  multiple  of  a  chi-square 
variate  with  2M  degrees  of  freedom.  The  probability  density  function  of  envelope 
Wj£  (reference  3,  equations  (3. 7-5)  and  (3. 7-10))  is  given  by 

Wk  (  "k\ 

P<wk>  -  —7  exP  ! - 2  •wk>0  ■  <4> 

<7  \  2«r  / 


Therefore  the  probability  density  function  of  z  (reference  4,  equations  (2.  l), 
(2. 3),  and  (2. 13))  is 


p(z)  = 


(M-l) 


/  VM-1  /  v 

hrfc)  axp  fi) 


,  z  >  0 


(5) 


The  probability  that  RV  z  exceeds  A  is 


Prob(z  > 


ro 


dz  p(z) 


*  tu-le* 
(M-l)! 


(6) 


the  last  step  via  repeated  integration  by  parts. 


The  probability  of  false  alarm,  PF’  and  the  probability  of  detection,  P^, 
are  obtained  by  using  (1)  in  (6): 


PF  =  exp(-A) 


M-  1 


£ 

m  =  0 


t 


(7) 

(8) 
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where  we  have  defined  normalized  threshold,  A,  as 


Again,  S/N  is  the  ratio  of  signal  power  to  noise  power  at  the  NBF  output  in 
figure  1  in  the  B  Hz  bandwidth;  it  is  not  the  spectral  level  in  a  1-Hz  band. 

Equation  (7)  can  be  solved  for  normalized  threshold  A  for  a  specified  PF 
and  number  of  independent  SE  samples  M.  (A  table  for  M  =  1(1)10,  16,  32, 
and  64,  and  PF  =  10"n,  n  =  1(1)8,  is  presented  in  Appendix  B. )  Then,  for  a 
specified  PD’  (8)  can  be  solved  for  the  required  values  of  S/N. 

RESULTS 

Curves  of  required  S/N  in  decibels  are  presented  in  figures  2  through  5 
for  PD  =  .  5,  .7,  .9,  and  .  99,  respectively.  Values  of  PF  equal  to  10“n, 
n  =  1(1)8,  are  considered  for  a  range  of  M  from  1  to  100.  (Since  the  curves 
were  evaluated  only  at  the  integer  points,  straight-line  interpolation  was  used 
for  ease  of  reading. )  M  is  to  be  interpreted  as  BT  +  1  in  figure  1  and  is  the 
equivalent  number  of  independent  SE  samples  in  observation  time  T, 

It  will  be  observed  from  figures  2  through  5  that  the  required  S/N  increases 
rapidly  as  M  decreases.  In  fact,  the  required  value  of  S/N  for  M  =  1  is 
given  by 

in  PF 

S/N=  T-—  -  !,(M=  !);  (10) 

in  PD 

this  result  follows  from  (7)  and  (8).  Thus  for  values  of  P^  near  unity,  a  very 
large  value  of  S/N  is  required  for  M  =  1. 

On  the  other  hand,  the  curves  are  fairly  similar  for  large  M,  there  being, 
for  example,  1.  9  dB  more  SNR  required  (at  M  -  100,  PF  -  10“^)  for  P^  =  .  9 
than  for  PD  =  .  5. 

The  difference  between  the  curves  decreases  as  PD  is  increased  in 
figures  2  through  5.  Part  of  this  difference  is  the  result  of  the  change  of  the 
ordinate  ncale  in  each  figure,  but  part  of  it  is  a  real  effect.  For  example,  the 
difference  between  the  curves  at  M  -  100,  PD  -  .  5  is  2. 1  dB,  whereas  the  dif¬ 
ference  between  the  curves  at  M  =  100,  Pp  -  .99  is  1.4  dB.  The  corresponding 
differences  at  M  =  1  are  3.  25  and  3.  0  dB,  respectively. 
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Figure  3.  Required  Signal-to-Noiee  Ratio  for  PD  =  .  7 
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It  is  anticipated  that  for  large  M  the  BV  z  in  (3)  could  be  approximated 
by  a  Gaussian  BV,  This  is  Investigated  in  Appendix  C,  where  the  approximate 
S/N  required  for  large  M  is  given  by  (C-8)  as 


+_1<PD) 

S/N  w - 

VS?  -  +"1  (PD) 


(11) 


where  4>_1(  )  is  the  Inverse  *(  )  function  (see  (C-5)>.  Equation  (11)  is  ex¬ 
pected  to  be  a  good  approximation  only  for  M  »  1,  and  only  for  Pq  signifi¬ 
cantly  less  than  a  point  that  is  discussed  further  in  Appendix  C. 


A  special  case  of  (11)  is  provided  by  PD  »  .  5  (see  (C-10)).  Converted  to 
decibels,  (11)  becomes 

10  log  (S/N)  «  10  log  (-  +“1  (PF))  -  5  log  M;  PD  *  .  5  ,  (12) 


which  has  the  familiar  S  log  M  (35  log  (BT))  decay  associated  with  energy 
detection  for  large  M,  Equation  (12)  is  plotted  in  figure  6  as  dashed  lines  for 
Pp  =  IQ"1,  10“2,  10-4,  and  10”8,  and  the  exact  results  are  plotted  as  solid 
lines.  The  values  of  the  additive  constant  in  (12)  are  given  in  table  1  (see  table 
26. 6  of  reference  5). 


Table  1.  Additive  Constant  in  (12) 


Pp 

10  log  (-+-1  (Pp)) 

10‘* 

1.08 

10“2 

3.67 

co 
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5.  70 

io-s 
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7. 16 

00 
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Figure  8.  Comparison  of  Exact  Results,  Gaussian  Approximation, 

and  Sine- Wave  Signal 
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Figure  6  shows  that  the  exact  results  approach  the  Gaussian  approximation 
asymptotically  as  M  increases.  However,  even  at  M=  100,  the  Gaussian 
approximation  is  optimistic  by  0.78  dB  for  Pp  =  10"8,  0.  53  dB  for  Pp  =  10"*, 
and  by  onjy  0. 19  dB  for  Pp  =  10”1. 

For  smaller  M,  such  as  M  =  10,  the  Gaussian  results  are  optimistic  by 
0.63,  1.07,  1.63,  and  2. 30  dB  at  PF  =  10-1,  10“2,  10”4,  and  10“8,  respec¬ 
tively,  but  for  M  *  1,  they  are  optimistic  by  2. 58,  3.  85,  5. 20,  and  6. 59  dB, 
respectively.  Thus  the  Gaussian  approximation,  and  (12)  in  particular,  should 
not  be  used  unless  M  is  very  large  compared  with  unity. 

For  other  values  of  PD,  (11)  can  be  converted  into  decibels  and  compared 
with  the  exact  results.  A  sample  calculation  for  PD  =  0.9,  Pp  =  10"®,  and 
M  =*  10  reveals  that  the  exact  S/N  required  is  6. 29  dB,  whereas  the  approxi¬ 
mation  (11)  indicates  5. 06  dB;  thus  the  Gaussian  approximation  is  optimistic  by 
1. 23  dB  in  this  case.  It  is  worth  noting,  however,  that  this  error  is  less  than 
the  corresponding  error  of  2. 00  dB  when  PD  =  0. 5,  Pp  =  10“6,  and  M  =  10. 

It  is  also  worth  noting  that  when  PD  *  0.  9,  Pp  =  10“8,  and  M  =  100,  the 
exact  S/N  Is  -1. 12  dB,  whereas  the  approximation  (11)  Indicates  -1.  60  dB; 
the  discrepancy  here  is  only  0.48  dB,  which  is  again  less  than  the  0. 67  dB  error 
for  PD  =  0. 5. 

Generally,  for  P^  /  0. 5,  the  approximation  (11)  does  not  plot  as  a  straight 
line  as  in  figure  6.  Rather,  the  plot  of  (11)  curves  up  sharply  as  M  decreases 
and,  in  fact,  goes  to  ~  at  M  =  [<t>-1  (Pd)]2-  t^us  approximation  (11)  can 
overestimate  the  required  S/N  for  small  enough  M  and  PD  >  0.  5.  Also,  as 
M  increases  from  small  values,  (11)  crosses  the  exact  S/N  curve  and  reaches 
a  maximum  undershoot  before  again  approaching  the  exact  S/N  curve  for 
large  M.  The  net  effect  is  that  (11)  is  not  a  good  approximation  to  use  unless 
both  M  »  1  and  M  »  Co-1  (PD)]2- 

If  the  signal  were  not  a  Gaussian  process,  but  a  sine  wave,  we  could  use  the 
results  of  reference  6  (equations  (B-21)  through  (B-25)  and  (013))  by  identifying 
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Equation  (13)  can  be  solved  for  S/N  In  terms  of  dp  (Figures  2  through  14 
in  reference  6  provide  the  desired  values  of  dp )  We  have  plotted  as  Xs  in 
figure  6  the  required  S/N,  in  decibels,  for  the  sine-wave  process  for  Py  * 
10"2  and  10-8.  Little  difference  is  evident  between  the  S/N  required  for  a 
Gaussian  signal  process  and  that  required  for  a  sine-wave-signal  process,  ex¬ 
cept  for  very  small  values  of  M.  The  sine-wave-signal  process  requires  a 
slightly  lower  value  of  S/N  than  does  the  Gaussian  signal  process  —  this 
conclusion  is  drawn  only  for  the  PD  =*  0. 5  case  considered  here. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Rules  of  thumb  for  determining  required  SNR,  such  as  those  developed  from 
a  Gaussian  approximation,  are  grossly  optimistic  predictions  when  the  number 
of  SE  samples,  M(=BT  +  1),  is  small.  Because  errors  of  many  decibels  can 
occur  for  M  of  the  order  of  unity,  the  results  shown  in  figures  2  through  5 
should  be  preferred  for  performance  predictions  in  narrowband  energy  detec¬ 
tion.  However,  these  results  should  still  be  considered  as  optimistic  because, 
in  practice,  unavoidable  center  frequency  and  bandwidth  mismatches  of  the 
analysis  filter  with  respect  to  the  signal  spectrum  occur  and  require  larger 
values  of  SNR  than  are  indicated  in  this  report.  (Bandwidth  mismatch  is  con¬ 
sidered  quantitatively  in  reference  1. ) 

The  false-alarm  and  detection  probabilities  themselves,  not  deflection 
criteria,  are  used  as  performance  measures  here.  Also,  no  assumptions  about 
the  absolute  or  relative  sizes  of  analysis  bandwidth  B  and  observation  time  T 
are  required;  only  their  product  is  important.  The  major  approximation  has 
been  the  replacement  of  BT  by  an  effective  number  of  independent  SE  samples, 
M-l.  This  approximation  is  difficult  to  be  precise  about,  without  an  exact 
analysis  as  conducted  in  reference  7.  (Details  of  the  exact  analysis  for  con¬ 
tinuous  detection  are  extremely  tedious,  as  Appendix  B  in  reference  7  attests. ) 
It  appears  that  the  next  order  of  business  should  be  an  exact  analysis  of  the 
sampled  system  in  figure  1  using  the  methods  in  references  8  and  9. 

The  results  in  this  report  are  actually  exact  for  a  rectangular  NBF,  if  BT 
is  an  integer;  then,  M  Is  precisely  BT  +  I.  This  follows  because  every  SE 
sample  in  figure  1  is  statistically  Independent  in  the  case  of  a  rectangu¬ 

lar  NBF,  and  there  are  T/B"1  +  1  samples  in  T  (see  Appendix  A). 
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Appendix  A 

EQUIVALENT  NUMBER  OF  INDEPENDENT  SQUARED-ENVELOPE  SAMPLES 


The  problem  here  is  to  find  on  equi  valent  number  of  independent,  squared- 
envelope  (SE)  samples  for  the  detector  in  figure  1  when  the  NBF  and  sampling 
increment  are  arbitrary  and  the  SE  samples  are  statistically  dependent.  With 
this  information,  the  results  of  this  report  can  be  extended  approximately  to 
other  signal  spectra,  filters,  and  sampling  plans. 


The  observation  interval  available  to  the  sampler  in  figure  1  is  (0,  T), 
during  which  time  an  arbitrary  number,  N,  of  equi-spaced  SE  samples  are 
taken.  The  time  between  SE  samples  is  then  At  =  T/(N- 1).  A  related  prob¬ 
lem  that  we  will  address  is:  What  value  of  At  should  be  chosen  to  make  the 
best  use  of  the  available  data? 

The  sum  of  N  (dependent)  SE  samples  in  figure  1  is  proportional  to 


u  * 


1 

2N 


N  - 1 


r 

k=0 


N-l 

S  frAt)  =  2N  L 
k  =  0 


(A-l) 


Using  the  probability  density  function  for  the  envelope  given  in  (4),  we  see  that 
the  mean  of  RV  u  is 


E{uf  =  i.E{e2(t)f=*2  .  (A-2) 

Now 

S2(t)=y^(t)  +  y^(t)  ,  (A-3) 

where  yc(t)  and  y8(t)  are  the  in-phase  and  quadrature  components  of  the  NBF 
output  y(t)  (see  reference  7,  Appendix  A).  If,  and  only  if,  the  spectrum  Gy (f) 
of  Gaussian  process  y(t)  is  symmetric  about  its  center  frequency  fQ,  then 
yc(t)  and  y8(t)  are  Independent  Gaussian  processes,  each  with  correlation 
function 


df  cos(2rfr)  G0(f),  Ro(0)  =  <r2 


(A-4) 


I 
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/ 

/ 


a 


\ 


i 

i 


where 


2Oy(f+fo).<>-fo| 

0  '  0,  f<-f„ 


(A-B) 


Then  each  SE  sample  is  the  sum  of  two  independent  samples. 
The  variance  of  u  follows  as 

v«)«|  -  k  S_M  (l  -  -I1)  Eo  (iTt) 


(A-6) 


Now  let  us  consider  a  new  BVv  formed  as  the  sum  of  M  statistically 
independent  SE  samples: 


1  “  2 
V=  2M  JE  wk 

k  =  l 


(A-7) 


RV  v  is  a  multiple  of  a  chi-square  variate  with  2M  degrees  of  freedom;  it  is 
also  a  scaled  version  of  (3).  We  wish  to  approximate  the  statistics  of  the  general 
RV  u  by  those  of  RV  v.  In  particular,  we  will  set  the  means  and  variances  of 
u  and  v  equal.  Using  (4)  and  (A-4),  we  find 


E{v[  =<r2, 

Varjvf=-i,*=iB>  •  (A-8) 


The  means  of  u  and  v  are  already  equal  and  the  variances  can  be  made  equal 
if  we  choose  M  in  (A-7)  (by  equating  (A-6)  and  (A-8))  according  to 


Equation  (A-9)  provides  a  definition  of  the  equivalent  number  of  independent  SE 
samples  in  RV  u.  (Since  the  right-hand  side  of  (A-9)  is  not  necessarily  an 
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integer,  M  could  be  chosen  as  the  closest  integer;  then  the  variances  of  u 
and  v  would  be  approximately  equal. )  In  the  special  case  when  At  is  chosen 
such  that  the  samples  of  g2(t)  in  (A-l)  are  all  uncorrelated,  then  Rq  (kAt)  *  0 
for  k  ji  0,  and  M  *  N. 

Two  special  cases  of  (A- 9)  are  worth  noting.  The  first  is  that  in  the  limit 
of  continuous  processing,  N— and  (A-9)  approaches 

1 

Mc  ■ - * - .  (A-10) 

/I  R  (Tx) 

dx(l-|x|)-~ - 

-1  R>) 


The  second  is  that  if  observation  time  T  is  large  compared  with  the  effective 
correlation  time  of  RQ(r),  then  (A-10)  becomes, 


MC=S 


r  .  ]2 

?  1 

/  *<>0(f) 

[ 

IV—  J  _  T 

ro  J 

£*«>  [•*> 

--  R  (0) 
o' 


where  we  have  used  (A-4),  (A-5),  and  (2),  and  where  B  is  the  effective  (or 
statistical)  bandwidth  of  the  positive-frequency  component  of  the  NBF  in  figure  1. 


Returning  to  the  general  result  for  the  equivalent  number  of  independent 
SE  samples,  (A-9),  we  note  from  figures  2  through  5  that  M  is  desired  to  be 
as  large  as  possible.  Therefore,  we  investigate  the  behavior  of  M  in  (A-9) 
as  a  function  of  N,  and  select  that  value  of  N  corresponding  to  a  maximum 
of  M;  the  best  N  will  not  be  infinite  (see  reference  10). 

We  must  investigate  (A-9)  through  particular  examples.  The  first  example 
is  a  rectangular  spectrum  for  the  NBF  output  y(t): 
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GyW  = 


™  *  lf±f0l<T^ 


G0(f)  =  {B  '» 

0,  otherwise 


Bc(r)  =  °2  SlTn^BrU  <r2  sine  (Br)  . 


Substitution  of  (A-12)  in  (A-9)  yields 


M  =■ 


(A-12) 


(A- 13) 


The  quantity  M  is  a  function  only  oftheBT  product  and  the  particular  value  of  N 
selected  and  is  plotted  versus  N  for  several  BT  products  in  figure  A-l.  It  will 
be  observed  that  M  increases  nearly  linearly  with  N  until  M  saturates  at  approx¬ 
imately  the  value  BT  +  1  when  N  =  BT  + 1.  The  corresponding  At  is  T(N  -  1)  « 
l/B;  that  ist  the  SE  samples  are  uncorrelated*  (RQ(k  At)  -  0  for  k  /  0)  when 
taken  at  the  bandwidth  rate  of  theNBF.  (Although  M  can  be  made  slightly  larger 
than  BT  +  1  by  choosing  N  one  larger  than  BT  +  1,  i.  e. ,  At  <  l/B,  the  gain  is 
negligible. )  Thus  for  the  rectangular  spectrum,  we  have  At  =  l/B,  M  =  BT  + 1. 


The  second  example  is  a  Gaussian  spectrum, 

G0(f)  =  *  exp(-2rf2/B2) 

R0(r)  =  a2exp(-  fB2r2)  . 


(A- 14) 


The  plot  of  M  versus  N  is  given  in  figure  A-2.  As  in  the  case  of  a  rectangu¬ 
lar  spectrum,  M  saturates  at  approximately  the  value  BT  +  1  when  N  equals 
BT  +  1,  and  At  =  l/B  is  again  the  recommended  sampling  interval.  However, 
the  SE  samples  are  no  longer  statistically  independent;  instead,  adjacent  SE 
samples  have  a  correlation  coefficient*  of  exp(-x)  =  0,043. 

♦The  correlation  coefficient  of  82(t)  is  R2  (r)/R^  (0). 
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NUMBER  OF  SQUARE  ^ENVELOPES  SAMPLES,  N 

Figure  A-l.  Equivalent  Number  of  Independent  Squared-Envelope 
Samples  for  Rectangular  Spectrum 


0  4  8  12  16  20  24  28 

NUMBER  OF  SQUARED-ENVELOPES  SAMPLES,  N 

Figure  A-2.  Equivalent  Number  of  Independent  Squared- Envelope 
Samples  for  Gaussian  Spectrum 
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The  third  example  is  a  Cauchy  spectrum. 


G0(f)  =  2,2 


(2  xf)2  +  B2 


(A-15) 


Rft(r)  =  *  exp(-B|r|)  . 


The  plot  of  M  versus  N  is  given  in  figure  A-3.  Although  M  again  saturates 
for  large  N  approximately  at  the  value  BT  +  1,  the  value  of  N  required  is 
larger  than  BT  +  1.  In  fact,  a  value  of  N  approximately  equal  to  2BT  is  re¬ 
quired  to  realize  the  saturation  value  of  M.  In  this  case,  the  sampling  increment 
At  equals  1/$B)  and  the  SE  samples  are  highly  correlated,  the  correlation 
coefficient  of  adjacent  SE  samples  being  1/e  =  0.  37. 


Thus,  for  three  widely  different  spectra,  the  equivalent  number  of  inde¬ 
pendent  SE  samples  is  given  approximately  by  BT  +  1,  where  B  is  the 
effective  bandwidth  of  the  NBF  output  spectrum.  The  sampling  increment 
should  be  approximately  l/B  for  fairly- rectangular  passband  characteristics, 
but  must  be  smaller  than  this  amount  for  rounded  filter  characteristics. 
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Appendix  B 
THRESHOLD  VALUES 

The  false  alarm  probability  is  given  in  (7).  Values  of  the  normalized  thresh¬ 
old  A  »  A  /{ShO  required  to  realize  specified  values  of  PF  are  given  in  table  B-l. 

Table  B-l.  Normalized  Threshold  Values  A 


M 

PF 

H 

1 

O 

H 

10“2 

10“3 

10"4 

io-6 

H* 

© 

1 

a 

10"7 

ESI 

B 

2. 303 

4.605 

6.908 

9.  210 

11.513 

13.  816 

16. 118 

18.421 

3.890 

6. 638 

9.233 

11.  756 

14.  237 

16.  688 

19. 120 

21. 536 

B 

5. 322 

8.406 

11.  229 

13.  928 

16. 554 

19. 129 

21. 669 

24.181 

B 

6.681 

10. 045 

13. 062 

15.  914 

18. 666 

21.  350 

23. 986 

26. 585 

5 

7.994 

11.  605 

14.  794 

17.  782 

20.648 

23.431 

26. 155 

28. 832 

6 

9.276 

13. 108 

16.455 

19.  567 

22.  538 

25. 413 

28.217 

30. 967 

7 

10.  532 

14. 571 

18. 062 

21.  290 

24.358 

27.  318 

30. 198 

33.016 

8 

11.  771 

16.  000 

19. 626 

22.  962 

26. 122 

29. 162 

32. 114 

34.997 

9 

12. 995 

17.403 

21. 156 

24. 595 

27.  841 

30.  957 

33. 977 

36. 922 

10 

14. 206 

18.  783 

22. 657 

26. 193 

29.  522 

32.  710 

35. 795 

38. 799 

16 

21.  292 

26. 743 

31.244 

35. 286 

30.047 

42.  616 

46.042 

49. 359 

32 

39. 430 

46. 608 

52. 358 

57. 417 

62. 052 

66.  394 

70. 519 

74.477 

64 

74. 443 

84. 067 

91.  593 

98. 106 

103. 996 

109.453 

114. 592 

119.482 

If  the  mean  value  of  decision  variable  z  in  figure  1  is  subtracted  before 
threshold  comparison,  the  number  of  standard  deviations  that  the  new  thresh¬ 
old  must  be  set  at  is  given  (for  signal  absent)  by 

A-  E{z[  A-M2N  A 


St.  Dev.  |z[  \/5T  2N  >/5T 


-  Vm1  . 


(B-l) 


The  results  in  table  B-l  may  be  substituted  in  equation  (B-l)  to  find  the  new 
thresholds.  For  example,  M  *  64,  Pp  =  10~3  requires  a  threshold  set  at 
3.449  standard  deviations. 
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Appendix  C 

GAUSSIAN  APPROXIMATION 


that 


For  large  M,  RVz  In  (3)  tends  toward  a  Gaussian  RV.  It  follows  from  (5) 


In  particular,, 


Then, 


*n«k(-  ^4-f,^  <*«*>“ 


B|z|  »M2»2, 

2 

Var  jz}  »  M{2*2)  . 


P<*)« 


V57ST  2<r2 


exp 


(a  -  M2<r2) 

- - g-,  , 

2M(2<r2) 


(C-l) 


(C-2) 


(C-3) 


and 


’vhere 


Prob(z  >  A)  w  ♦ 


(vh,-^F  p)  ’ 


(C-4) 


\ 


A 

♦  (x)  •  J  dt  (2x)'l/2  exp(-t2/2)  . 


(C-5) 


An  immediate  problem  of  the  Gaussian  approximation  is  indicated  by  (C-4): 
as  <r2  —  the  right-hand  side  of  (C-4)  tends  to  ♦(VST),  which  is  less  than  1. 
Thus  arbitrarily  high  detection  probabilities  can  never  be  realized  by  die 
Gaussian  approximation  unless  M  »  1.  This  is  not  true  of  the  exact  distribu¬ 
tion  of  z  for  any  M,  as  (6)  ahowa  that  the  Prob(z  >  A)  —  1  as  <r2  — 

The  approximations  to  Pp  and  PD  are  obtained  by  using  (1)  and  (9)  in 
(C-4): 
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(C-6) 

(C-7) 


Equations  (C-6)  and  (C-7)  may  be  solved  for  S/N  (by  eliminating  A)  to  obtain 


+"1  (PD)  -  +’1  (PF) 
V®-*  1  (PD) 


(C-8) 


where  «l>“1(  )  is  the  inverse  <t>(  )  function.  (The  impossibility  of  requiring 
PD  >  <t>(VM)  is  now  made  obvious  by  the  denominator  of  (C-8). ) 

For  large  M,  (C-8)  decays  as  M”*^,  regardless  of  the  choice  of 
and  Pp.  However,  the  point  at  which  this  behavior  develops  depends  upon  how 
large  PD  is. 

A  special  case  of  (C-8)  is  provided  by  Pp  =  .  5;  then  *“1(PD)  =  0  and 

-*-1  (PF) 

S/N  « - - ;  PD  -  .  5  .  (C-9) 

Converted  to  decibels,  (C-9)  is 

10  log  (S/N) «  10  log  (-  *“1  (PF))  -  5  log  M;  PD  »  .  6  .  (C-10) 
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